Analisis de ciclos en (LAS BAMBAS)¶

Import Packages¶

Este chunk configura el entorno de análisis de datos y las librerías necesarias para el proyecto. Se establecen las configuraciones de visualización y se importan todas las dependencias requeridas para:

📊 Análisis de datos: pandas, numpy para manipulación y cálculos estadísticos 📈 Visualización: plotly para gráficos interactivos con configuración optimizada para notebooks 🧮 Estadísticas: scipy.stats, statsmodels, pingouin para pruebas estadísticas y análisis avanzados 🔬 Análisis avanzado: gaussian_kde para estimación de densidad kernel y cálculos de probabilidad

La configuración incluye supresión de warnings para mantener limpia la salida y ajustes de plotly para garantizar compatibilidad con exportación HTML estática.

Import data and preprocessing¶

Este chunk se encarga de la carga y transformación inicial de los datos de rendimiento TKPH (Toneladas por Kilómetro por Hora) de los camiones mineros. Se realiza un preprocesamiento completo que incluye:

🚛 Filtrado de modelos: Se seleccionan únicamente los camiones de interés:

  • KOM 980E (Komatsu 980E)
  • CAT 797F (Caterpillar 797F)
  • CAT 798AC (Caterpillar 798AC)

⚙️ Creación de variables calculadas:

  • Velocidad: Cálculo de velocidad promedio del ciclo (km/h)
  • Hora: Extracción de la hora del día (0-23) para análisis temporal
  • Peso delantero ajustado: Cálculo del peso distribuido en las llantas delanteras según especificaciones técnicas de cada modelo

🎯 Reglas de negocio aplicadas:

  • Porcentajes de distribución de peso específicos por modelo de camión
  • Ajuste especial del +6% para operaciones desde pit CHALCOBAMBA hacia Chancador
  • Manejo de valores faltantes y validación de datos

Los datos procesados quedan listos para análisis estadísticos, visualizaciones y generación de reportes de performance operacional.

Especificaciones Técnicas - Distribución de Peso en Llantas Delanteras

MODELO EVW_FRONT GVW_FRONT Descripción
CAT 797F 48.5% 32.8% Caterpillar 797F
CAT 798AC 47.0% 33.0% Caterpillar 798AC
KOM 980E 48.5% 32.8% Komatsu 980E
📏 1. Calculando Velocidad...

📏 2. Calculando Hora...

📏 3. Peso ajustado...

Analisis de payload¶

Este chunk realiza un análisis estadístico completo de la distribución de tonelaje de carga (payload) utilizando dos enfoques complementarios:

Análisis de Distribución por Rangos Operacionales¶

Se implementa una clasificación por rangos basada en umbrales operacionales críticos:

  • Verde (< 363 ton): Zona de subcarga
  • Amarilla (363-400 ton): Zona de transición
  • Roja (400-436 ton): Zona objetivo
  • Negra (> 436 ton): Zona de sobrecarga

El análisis genera un gráfico dual que combina:

  1. Histograma de densidad: Visualización de la distribución real de tonelaje
  2. Tabla resumen: Porcentaje de ciclos que caen en cada rango operacional

Cálculo de Probabilidades usando Kernel Density Estimation (KDE)¶

Se calcula P(Tonelaje > 400) mediante dos métodos estadísticos:

Método KDE: Utiliza estimación de densidad kernel para crear una función de densidad continua, seguida de integración numérica (trapezoides y cuadratura adaptativa de scipy) para calcular probabilidades exactas.

Método Empírico: Cálculo directo basado en la proporción observada de datos que superan el umbral.

La comparación entre ambos métodos valida la robustez del análisis. El chunk incluye visualización interactiva que muestra el área bajo la curva correspondiente a la probabilidad calculada, proporcionando una interpretación intuitiva del resultado estadístico.

📊 CREANDO GRÁFICO DE DENSIDAD DE TONELAJE
==================================================
✅ Datos de tonelaje: 68061 registros válidos
Rango: 321.0 - 464.0 toneladas

📋 DISTRIBUCIÓN POR RANGOS
==================================================
📈 ESTADÍSTICAS ADICIONALES
==================================================
Total de ciclos analizados: 68,061
Tonelaje promedio: 400.4 ton
Tonelaje mediano: 401.0 ton
Desviación estándar: 15.8 ton
🔬 CÁLCULO DE P(TONELAJE > 400) USANDO KERNEL DENSITY ESTIMATION
======================================================================
📊 Datos: 68,061 observaciones de tonelaje
🎯 Calculando: P(Tonelaje > 400)

📈 PASO 1: CREANDO KERNEL DENSITY ESTIMATOR
==================================================
✅ KDE creado exitosamente
Bandwidth (ancho de banda): 0.1080
Número de puntos de datos: 68061
Rango de datos: [321.0, 464.0] toneladas

🧮 MÉTODO 1: INTEGRACIÓN POR TRAPEZOIDES
=============================================
Puntos de integración: 2000
Rango de integración: [400, 464.0]
✅ P(Tonelaje > 400) = 0.514271
   En porcentaje: 51.427%

🔬 MÉTODO 2: INTEGRACIÓN SCIPY (CUADRATURA ADAPTATIVA)
============================================================
✅ P(Tonelaje > 400) = 0.514271
   En porcentaje: 51.427%
   Error estimado: 1.08e-12

📊 COMPARACIÓN CON MÉTODO EMPÍRICO
=============================================
Empírico:     P(T > 400) = 0.501095 (50.109%)
KDE (Scipy quad): P(T > 400) = 0.514271 (51.427%)

Diferencia absoluta: 0.013177
Diferencia porcentual: 1.318 puntos porcentuales
✅ Buena concordancia entre métodos

Ciclos observados > 400: 34105 de 68061 total

📊 CREANDO VISUALIZACIÓN...
💻 FUNCIÓN REUTILIZABLE PARA KDE
========================================

🔧 EJEMPLO DE USO:
--------------------
Método: KDE + Scipy quad
P(Tonelaje > 400) = 0.514271
Porcentaje: 51.427%
Error estimado: 1.08e-12

✅ RESULTADOS FINALES
==============================
🎯 Umbral analizado: 400 toneladas
📊 Método usado: KDE con Scipy quad
🔢 P(Tonelaje > 400) = 0.514271
📈 Porcentaje: 51.427%
📋 Interpretación: De cada 1000 ciclos, aproximadamente 514 tendrán tonelaje > 400

Pesos delanteros¶

Este chunk implementa un análisis exhaustivo de la regla operacional especial que aplica un ajuste del +6% al peso distribuido en llantas delanteras para operaciones desde el pit CHALCOBAMBA hacia destinos que contienen "Chancador".

Análisis de Impacto Operacional¶

Se cuantifica la frecuencia de aplicación de la regla especial:

  • Conteo total de operaciones afectadas
  • Porcentaje de ciclos con ajuste vs operaciones normales
  • Análisis de rutas específicas desde CHALCOBAMBA

Visualización de Distribuciones Comparativas¶

Genera un dashboard de 4 paneles que compara las distribuciones de peso delantero:

  1. Histogramas de densidad: EVW y GVW delantero separados por condición normal/ajustada
  2. Box plots comparativos: Visualización de estadísticas descriptivas y detección de outliers

Análisis Estadístico Robusto¶

Implementa pruebas estadísticas no paramétricas (Mann-Whitney U) para determinar si las diferencias entre grupos son estadísticamente significativas, evitando asumir normalidad de los datos.

Segmentación Multi-dimensional¶

Desglosa el análisis por:

  • MODELO: Distribución de ajustes por tipo de camión (KOM 980E, CAT 797F, CAT 798AC)
  • EQUIPO: Identifica equipos específicos más afectados por la regla, ordenados por número absoluto de ajustes

Visualizaciones Multi-panel¶

  • Gráficos de barras apiladas mostrando proporción de ajustes por modelo
  • Box plots comparativos de peso delantero por modelo
  • Heatmap por equipo correlacionando ciclos totales con porcentaje de ajustes
🎯 ANÁLISIS COMPLETO - REGLA ESPECIAL AJUSTE CHALCOBAMBA
======================================================================
✅ Todas las columnas necesarias están presentes

📊 ESTADÍSTICAS GENERALES 
=======================================================
Tipo_Registro Cantidad Porcentaje (%) Media_Distancia_KM
0 Con Ajuste 16,485 24.22% 16.62
1 Sin Ajuste 51,576 75.78% 10.29
2 Total 68,061 100.00% 11.83
🏗️ ANÁLISIS POR PIT Y DESTINO
================================================================================
--------------------------------------------------------------------------------

📊 RESUMEN COMPLETO POR PIT:
------------------------------------------------------------------------------------------------------------------------
Num_Ciclos Dist_Media_KM Dist_Total_KM Dist_Min_KM Dist_Max_KM Dist_Std_KM Ciclos_Con_Ajuste Horas_Totales Pct_Con_Ajuste Pct_Del_Total
PIT
CHALCOBAMBA 34240 11.24 384863.48 0.29 61.34 6.64 16485 23575.57 48.1 50.31
FERROBAMBA 33821 12.42 420065.65 0.26 55.90 5.21 0 29689.24 0.0 49.69
📋 DESTINOS DESDE CHALCOBAMBA:
Num_Ciclos Dist_Media Dist_Min Dist_Max Dist_Std Con_Ajuste Pct_Con_Ajuste
COMP_DESTINO
Botadero 22,755 8.18 0.54 43.51 4.50 6,967 30.6%
Chancador 5,055 21.05 0.29 61.34 5.40 5,055 100.0%
Stock Sulfuros 4,685 14.27 0.88 46.73 4.41 3,349 71.5%
Stock Mixtos 1,525 15.53 2.22 53.28 4.78 985 64.6%
Inpit 208 7.58 0.49 22.69 4.15 124 59.6%
Stock Óxidos 12 12.30 7.98 15.78 2.85 5 41.7%
📈 CREANDO GRÁFICOS DE DENSIDAD DE PESO DELANTERO...
🧮 ANÁLISIS ESTADÍSTICO COMPARATIVO
=============================================

PESO GVW DELANTERO:
-------------------------
Normal    (n=51,576): μ=110.5, σ=2.8
Ajustado  (n=16,485): μ=114.6, σ=3.0
Diferencia de medias: +4.1 ton (+3.7%)

📊 PRUEBAS ESTADÍSTICAS:
-------------------------
GVW - Mann-Whitney U: p-value = 0.00e+00
GVW - Diferencia significativa: SÍ

🚛 ANÁLISIS POR MODELO DE CAMIÓN
=============================================
RESUMEN POR MODELO:
-------------------------
MODELO Total_Ciclos Con_Ajuste Sin_Ajuste Pct_Ajuste GVW_Delantero_Prom
1 CAT 797F 33619 7914 25705 23.54 111.89
0 KOM 980E 32750 8274 24476 25.26 111.24
2 CAT 798AC 1692 297 1395 17.55 109.20
🔧 ANÁLISIS POR EQUIPO
==============================
TOP 10 EQUIPOS POR NÚMERO DE CICLOS AJUSTADOS:
----------------------------------------
HT112           | KOM 980E   | 2,462 ciclos |  653 ajustes ( 26.5%)
HT113           | KOM 980E   | 2,436 ciclos |  635 ajustes ( 26.1%)
HT083           | CAT 797F   | 2,410 ciclos |  625 ajustes ( 25.9%)
HT080           | CAT 797F   | 2,459 ciclos |  625 ajustes ( 25.4%)
HT104           | KOM 980E   | 2,263 ciclos |  620 ajustes ( 27.4%)
HT108           | KOM 980E   | 2,307 ciclos |  615 ajustes ( 26.7%)
HT115           | KOM 980E   | 2,276 ciclos |  610 ajustes ( 26.8%)
HT101           | KOM 980E   | 2,398 ciclos |  599 ajustes ( 25.0%)
HT071           | CAT 797F   | 2,304 ciclos |  588 ajustes ( 25.5%)
HT073           | CAT 797F   | 2,410 ciclos |  579 ajustes ( 24.0%)
✅ RESUMEN EJECUTIVO - REGLA AJUSTE CHALCOBAMBA
============================================================
📊 IMPACTO OPERACIONAL:
• Total de operaciones analizadas: 68,061
• Operaciones con ajuste especial: 16,485 (24.22%)

⚖️ IMPACTO EN PESO DELANTERO:
• GVW Delantero: +4.1 ton promedio (+3.7%)

🚛 DISTRIBUCIÓN POR MODELO:
• CAT 797F: 7914 ajustes de 33,619 total (23.5%)
• KOM 980E: 8274 ajustes de 32,750 total (25.3%)
• CAT 798AC: 297 ajustes de 1,692 total (17.6%)

📉 ANÁLISIS DE PESOS TRASEROS
⏱️ ANÁLISIS FULL_TRAVEL (HORAS)
Promedio FULL_TRAVEL_HRS Normal:   0.38 hrs
Promedio FULL_TRAVEL_HRS Ajustado: 0.39 hrs
⏱️ ANÁLISIS POR TIEMPO - MODELO Y AJUSTE
============================================================

📊 RESUMEN POR MODELO Y TIPO DE AJUSTE:
------------------------------------------------------------------------------------------
MODELO Tipo_Ajuste Num_Ciclos Horas_Totales Horas_Promedio_Ciclo Dist_Media_KM Pct_Tiempo_Total
0 CAT 797F Sin Ajuste 25705 18739.52 0.73 10.15 35.18
1 CAT 797F Con Ajuste (+6%) 7914 7334.56 0.93 16.52 13.77
2 CAT 798AC Sin Ajuste 1395 1047.91 0.75 11.06 1.97
3 CAT 798AC Con Ajuste (+6%) 297 252.96 0.85 15.97 0.47
4 KOM 980E Sin Ajuste 24476 18367.39 0.75 10.40 34.48
5 KOM 980E Con Ajuste (+6%) 8274 7522.47 0.91 16.74 14.12
📈 COMPARACIÓN DE TIEMPOS PROMEDIO:
------------------------------------------------------------

CAT 797F:
  Sin ajuste:  0.73 hrs/ciclo
  Con ajuste:  0.93 hrs/ciclo
  Diferencia:  +0.20 hrs (+27.4%)

CAT 798AC:
  Sin ajuste:  0.75 hrs/ciclo
  Con ajuste:  0.85 hrs/ciclo
  Diferencia:  +0.10 hrs (+13.3%)

KOM 980E:
  Sin ajuste:  0.75 hrs/ciclo
  Con ajuste:  0.91 hrs/ciclo
  Diferencia:  +0.16 hrs (+21.3%)

============================================================
✅ Análisis temporal completo
Archivo guardado: analisis_tiempo_por_modelo.html
📊 MATRIZ PIT vs DESTINO (Número de operaciones):
----------------------------------------------------

🔍 IDENTIFICANDO DESTINOS COMPARTIDOS...
PITs encontrados: CHALCOBAMBA, FERROBAMBA

✅ RESULTADO FINAL:
====================================================================================================
DESTINO Num_PITs PITs_Str Total_Operaciones Distancia_Promedio Coef_Variacion_Dist pct_CHALCOBAMBA pct_FERROBAMBA
0 CHANCADO 2 FERROBAMBA | CHALCOBAMBA 11391 15.05 47.36 44.377 55.623
2 BOT-F02-4195-RAM 2 FERROBAMBA | CHALCOBAMBA 3489 15.39 33.48 0.717 99.283
4 BOT-4525-PD-HHROESTE 2 CHALCOBAMBA | FERROBAMBA 670 4.14 94.21 99.851 0.149
1 BOT-SUR-4170 2 FERROBAMBA | CHALCOBAMBA 177 15.16 34.20 1.130 98.870
3 BOT-ZARANDA-ING-RECLA 2 FERROBAMBA | CHALCOBAMBA 53 17.93 23.55 3.774 96.226
5 IP-PISO-P12 2 FERROBAMBA | CHALCOBAMBA 32 6.02 93.45 46.875 53.125
6 IP-PISO-LD001 2 CHALCOBAMBA | FERROBAMBA 4 6.56 84.60 50.000 50.000

Analisis temporal¶

⏰ Análisis por Hora del Día¶

  • Se analizaron los ciclos en función de la hora del día.
  • Se compararon:
    • Volumen de ciclos realizados
    • Tonelaje transportado (mediana por hora)
    • Tiempo de ciclo (mediana en horas)
    • Distancia recorrida (mediana en km)
  • Identificación de horas pico y valle de actividad.

🔍 Análisis de Gaps (huecos operativos)¶

  • Un gap se define como el tiempo transcurrido entre el fin de un ciclo (dumping) y el inicio del siguiente (loadloading).
  • Se midieron:
    • Cantidad de gaps por hora y por equipo
    • Duración total y promedio de los gaps
  • Visualizaciones clave:
    • Distribución global de gaps por hora del día.
    • Heatmap global de gaps (equipos vs horas).
    • Heatmap diario filtrable por fecha (permite revisar jornadas específicas).
    • Serie temporal del % de tiempo total perdido en gaps diarios.

🕒 Análisis de Downtime (tiempo muerto)¶

  • El tiempo muerto se calculó como la proporción de espera (colas, spot, descarga) sobre el tiempo total de ciclo.
  • Indicadores generados:
    • Downtime promedio por ciclo.
    • Downtime ponderado por duración de ciclo.
    • Horas totales de downtime acumuladas.
  • Visualizaciones:
    • Boxplot comparativo de downtime % por modelo de equipo.
    • Ranking de las 20 combinaciones FECHA–EQUIPO con mayor downtime ponderado.
⏰ ANÁLISIS DE OPERACIONES POR HORA DEL DÍA
==================================================
✅ Columna 'Hora' encontrada

📊 CREANDO ANÁLISIS HORARIO...
RESUMEN POR HORA:
--------------------------------------------------------------------------------
Hora |  Ciclos | Tonelaje | Tiempo(h) | Distancia(km)
--------------------------------------------------------------------------------
   0 |    3251 |    400.0 |      0.77 |         11.4
   1 |    3147 |    399.0 |      0.78 |         11.7
   2 |    2816 |    400.0 |      0.81 |         11.9
   3 |    2943 |    400.0 |      0.81 |         12.0
   4 |    3037 |    400.0 |      0.77 |         11.3
   5 |    2674 |    401.0 |      0.77 |         11.4
   6 |     890 |    399.0 |      0.68 |          9.6
   7 |    1131 |    403.0 |      0.76 |         11.0
   8 |    2656 |    401.0 |      0.76 |         11.8
   9 |    3032 |    402.0 |      0.78 |         12.0
  10 |    3260 |    401.0 |      0.74 |         11.5
  11 |    3301 |    401.0 |      0.74 |         11.3
  12 |    3163 |    401.0 |      0.78 |         11.8
  13 |    3015 |    400.0 |      0.81 |         12.2
  14 |    3131 |    401.0 |      0.80 |         12.1
  15 |    3233 |    402.0 |      0.74 |         11.0
  16 |    3477 |    401.0 |      0.76 |         11.4
  17 |    3461 |    401.0 |      0.73 |         10.7
  18 |    1279 |    398.0 |      0.60 |          8.3
  19 |    1874 |    401.0 |      0.70 |          9.8
  20 |    3392 |    400.0 |      0.74 |         11.1
  21 |    3315 |    400.0 |      0.75 |         11.5
  22 |    3309 |    400.0 |      0.76 |         11.6
  23 |    3274 |    400.0 |      0.74 |         11.1
📈 ESTADÍSTICAS HORARIAS ADICIONALES
==================================================
ACTIVIDAD POR HORAS:
• Hora pico: 16h con 3477 ciclos
• Hora valle: 6h con 890 ciclos

TONELAJE POR HORAS:
• Mayor tonelaje: 403.0 ton a las 7h
• Menor tonelaje: 398.0 ton a las 18h
• Variación: 5.0 ton (1.3%)
⏰ Análisis temporal por horas completado!
📅 ANÁLISIS DE OPERACIONES POR MES
==================================================
🔄 Creando columna 'Mes'...
✅ Columna 'Mes' creada desde FECHA
✅ Columna 'Mes' encontrada

📊 CREANDO ANÁLISIS MENSUAL...

RESUMEN POR MES:
------------------------------------------------------------------------------------------------------------------------
         Mes |   Ciclos |  Tonelaje |  Total Ton |  Tiempo(h) |  Distancia |  Ajustes | % Ajustes
------------------------------------------------------------------------------------------------------------------------
       Enero |     7825 |     398.0 |    3105000 |       0.74 |        9.9 |     1958 |      25.0%
     Febrero |     6404 |     398.0 |    2548231 |       0.82 |       10.3 |     1092 |      17.1%
       Marzo |     8372 |     401.0 |    3355426 |       0.74 |        9.4 |     1479 |      17.7%
       Abril |     7971 |     402.0 |    3203266 |       0.66 |        8.8 |     1549 |      19.4%
        Mayo |     8042 |     401.0 |    3225433 |       0.80 |       12.1 |     2182 |      27.1%
       Junio |     8116 |     402.0 |    3259431 |       0.80 |       12.8 |     2591 |      31.9%
       Julio |     7290 |     401.0 |    2924929 |       0.78 |       12.6 |     1847 |      25.3%
      Agosto |     8221 |     401.0 |    3296542 |       0.74 |       12.1 |     2122 |      25.8%
  Septiembre |     5820 |     401.0 |    2334622 |       0.75 |       12.5 |     1665 |      28.6%
📈 ESTADÍSTICAS MENSUALES ADICIONALES
==================================================
ACTIVIDAD POR MESES:
• Mes pico: Marzo con 8372 ciclos
• Mes valle: Septiembre con 5820 ciclos
• Variación: 2552 ciclos (43.8%)

TONELAJE POR MESES:
• Mayor tonelaje total: 3,355,426 ton en Marzo
• Menor tonelaje total: 2,334,622 ton en Septiembre
• Variación: 1,020,804 ton (43.7%)

AJUSTES POR MESES:
• Mayor % de ajustes: Junio con 31.9%
• Menor % de ajustes: Febrero con 17.1%

PROMEDIOS MENSUALES:
• Ciclos promedio por mes: 7562
• Tonelaje promedio por mes: 3,028,098 ton
• % promedio de ajustes: 24.2%

TENDENCIAS SEMESTRALES:
• Primer semestre - Ciclos promedio: 7788
• Segundo semestre - Ciclos promedio: 7110
• Tendencia anual: a la baja
📅 Análisis temporal por meses completado!
📈 ESTADÍSTICAS DE TENDENCIA
==================================================
• Pendiente de la tendencia: 0.0098% por día
• Coeficiente de correlación (R): 0.122
• Coeficiente de determinación (R²): 0.015
• P-valor: 0.047908
• Tendencia estadísticamente significativa: SÍ

📊 COMPARACIÓN TEMPORAL:
• Promedio primer tercio: 52.99%
• Promedio último tercio: 54.81%
• Cambio: +3.4%

💡 RESUMEN EJECUTIVO:
• Total de equipos analizados: 31
• Período analizado: 2025-01-01 a 2025-09-20
• Promedio diario de tiempo perdido: 54.53%
• Máximo tiempo perdido en un día: 67.01%
• Mínimo tiempo perdido en un día: 7.01%
📊 ESTADÍSTICAS DETALLADAS POR MODELO
======================================================================

🔧 KOM 980E:
   • Número de ciclos: 32,750
   • Media: 11.25%
   • Mediana: 8.63%
   • Desviación estándar: 9.26%
   • Q1 (25%): 4.35%
   • Q3 (75%): 15.50%
   • Min-Max: 0.09% - 68.68%
   • Ciclos con downtime > 50%: 121 (0.4%)

🔧 CAT 797F:
   • Número de ciclos: 33,619
   • Media: 11.77%
   • Mediana: 9.14%
   • Desviación estándar: 9.38%
   • Q1 (25%): 4.65%
   • Q3 (75%): 16.25%
   • Min-Max: 0.08% - 66.20%
   • Ciclos con downtime > 50%: 129 (0.4%)

🔧 CAT 798AC:
   • Número de ciclos: 1,692
   • Media: 11.99%
   • Mediana: 9.70%
   • Desviación estándar: 9.40%
   • Q1 (25%): 4.93%
   • Q3 (75%): 16.48%
   • Min-Max: 0.12% - 60.27%
   • Ciclos con downtime > 50%: 8 (0.5%)
🔍 ANÁLISIS COMPARATIVO ENTRE MODELOS
==================================================
• Modelo con MENOR downtime promedio: KOM 980E (11.2%)
• Modelo con MAYOR downtime promedio: CAT 798AC (12.0%)
• Diferencia entre modelos: 0.7 puntos porcentuales
• Modelo MÁS consistente: KOM 980E (Std: 9.3%)
• Modelo MENOS consistente: CAT 798AC (Std: 9.4%)

Analisis TKPH¶

ANÁLISIS TKPH COMPLETO - CON KDE 3D Y FILTRO DE HORAS
=================================================================
DATOS ANTES DEL FILTRADO DE HORAS
Total de registros: 7,658
Período: 2025-01-01 00:00:00 a 2025-09-19 00:00:00
Equipos únicos: 31
Rango de horas: 0.22 - 17.21

ANÁLISIS KDE 3D - DATOS ORIGINALES
========================================
Datos para KDE 3D: 7,658 registros
APLICANDO FILTRO: HORAS >= 1
===================================
IMPACTO DEL FILTRO:
• Registros antes: 7,658
• Registros con horas < 1: 149 (1.9%)
• Registros después: 7,509
• Registros eliminados: 149 (1.9%)

KDE 3D - DATOS FILTRADOS (Horas >= 1)
========================================
Datos para KDE 3D filtrado: 7,509 registros
ESTADÍSTICAS DESCRIPTIVAS - DATOS FILTRADOS
=======================================================
       TKPH_Original  TKPH_Operacional  Hours_operation
count        7509.00           7509.00          7509.00
mean         1158.46           1334.52             7.08
std           178.68            191.29             2.99
min           409.44            530.31             1.00
25%          1044.24           1205.28             4.92
50%          1159.97           1332.09             6.99
75%          1273.94           1460.43             9.17
max          2174.10           2364.78            17.21

DIFERENCIAS (DATOS FILTRADOS):
• Diferencia promedio: 176.06 TKPH
• Diferencia porcentual promedio: 15.6%
• Correlación entre métodos: 0.962
======================================================================
📊 ANÁLISIS: CORRELACIÓN TKPH vs PRECIPITACIÓN
======================================================================

COEFICIENTES DE CORRELACIÓN (Pearson):
• TKPH Original vs Precipitación: -0.884
• TKPH Operacional vs Precipitación: -0.882

INTERPRETACIÓN:
• TKPH Original: Correlación Muy fuerte
• TKPH Operacional: Correlación Muy fuerte

RESUMEN MENSUAL:
------------------------------------------------------------------------------------------
            TKPH Original (Prom)  TKPH Operacional (Prom)  Precipitación (mm)
Nombre_Mes                                                                   
April                    1074.32                  1240.79                  26
August                   1273.35                  1459.58                   1
February                 1024.46                  1177.87                  48
January                  1110.80                  1287.06                  25
July                     1199.50                  1395.33                   1
June                     1266.67                  1448.97                   1
March                    1009.82                  1177.98                  64
May                      1178.99                  1346.36                   3
September                1346.91                  1539.37                   2

ANÁLISIS DE TEMPORADA DE LLUVIAS:
--------------------------------------------------

TEMPORADA DE LLUVIAS (>30mm):
Meses: February, March
• TKPH Original promedio: 1017
• TKPH Operacional promedio: 1178
• Precipitación promedio: 56.0 mm

TEMPORADA SECA (≤10mm):
Meses: August, July, June, May, September
• TKPH Original promedio: 1253
• TKPH Operacional promedio: 1438
• Precipitación promedio: 1.6 mm

IMPACTO DE LA PRECIPITACIÓN:
• Diferencia TKPH Original (Seca - Lluvias): +236 (+23.2%)
• Diferencia TKPH Operacional (Seca - Lluvias): +260 (+22.1%)

======================================================================
✅ Gráfico guardado: tkph_precipitacion_mensual.html
======================================================================

PRUEBAS DE HIPÓTESIS ESTADÍSTICAS - DATOS FILTRADOS
============================================================
PRUEBAS DE NORMALIDAD (Shapiro-Wilk):
• TKPH Original: W = 0.9930, p-valor = 6.18e-15
• TKPH Operacional: W = 0.9973, p-valor = 1.06e-07
• Distribución normal: NO (α = 0.05)

PRUEBA DE DIFERENCIAS PAREADAS:
• Prueba utilizada: Wilcoxon signed-rank
• W-estadístico: 27341.0000
• P-valor: 0.00e+00
• Diferencia estadísticamente significativa: SÍ (α = 0.05)

TAMAÑO DEL EFECTO:
• Cohen's d: 3.333
• Magnitud del efecto: grande

ANÁLISIS MENSUAL DETALLADO - DATOS FILTRADOS
=======================================================
January:
  • N = 867, Diferencia promedio = 176.26
  • P-valor = 0.0000, Significativo: SÍ
February:
  • N = 800, Diferencia promedio = 153.40
  • P-valor = 0.0000, Significativo: SÍ
March:
  • N = 891, Diferencia promedio = 168.16
  • P-valor = 0.0000, Significativo: SÍ
April:
  • N = 862, Diferencia promedio = 166.47
  • P-valor = 0.0000, Significativo: SÍ
May:
  • N = 894, Diferencia promedio = 167.37
  • P-valor = 0.0000, Significativo: SÍ
June:
  • N = 867, Diferencia promedio = 182.30
  • P-valor = 0.0000, Significativo: SÍ
July:
  • N = 895, Diferencia promedio = 195.83
  • P-valor = 0.0000, Significativo: SÍ
August:
  • N = 890, Diferencia promedio = 186.23
  • P-valor = 0.0000, Significativo: SÍ
September:
  • N = 543, Diferencia promedio = 192.46
  • P-valor = 0.0000, Significativo: SÍ

ESTADÍSTICAS DEL ANÁLISIS KDE 3D
========================================
KDE 3D DATOS ORIGINALES:
• TKPH Original - Densidad máxima: 0.00028106
• TKPH Operacional - Densidad máxima: 0.00025747

KDE 3D DATOS FILTRADOS:
• TKPH Original - Densidad máxima: 0.00028728
• TKPH Operacional - Densidad máxima: 0.00026526

PUNTOS DE MÁXIMA DENSIDAD (DATOS FILTRADOS):
TKPH Original:
• TKPH: 1179
• Horas: 6.8
TKPH Operacional:
• TKPH: 1330
• Horas: 6.8

RESUMEN EJECUTIVO - ANÁLISIS COMPLETO
=============================================
IMPACTO DEL FILTRO DE HORAS:
• Se eliminaron 149 registros (1.9%)
• Análisis final con 7,509 registros

RESULTADOS FINALES (Horas >= 1):
• El TKPH Operacional es 176.1 unidades mayor que el TKPH Original
• Esta diferencia ES estadísticamente significativa
• Correlación entre métodos: 0.962
• Fuerza de correlación: Fuerte
• El tamaño del efecto es grande

INSIGHTS CLAVE:
• 1.9% de los datos originales tenían < 1 hora de operación
• TKPH promedio de registros eliminados: Orig=1137.8, Op=1289.7
• Concentración máxima de densidad KDE:
  - TKPH Original: 0.00028728
  - TKPH Operacional: 0.00026526
  - TKPH Original muestra mayor concentración de valores

📦 BOXPLOTS SOLO TKPH_Original - POR AJUSTE
=======================================================